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Abstract 

We propose a solution to the 'cuspy-core' problem by extending the geodesic equations of motion 
using the Dark Energy length scale Xde = c/(A£)£;G)^/^. This extension does not affect the motion 
of photons; gravitational lensing is unchanged. A cosmological check of the theory is made, and erg 
is calculated to be 0.68±o.ii) compared to 0.76llQ'Q4g for WMAP. We estimate the fractional density 
of matter that cannot be determined through gravity at 0.197±o.oi7) compared to 0.196I'Iq'q26, the 
fractional density of nonbaryonic matter. The fractional density of matter that can be determined 
through gravity is estimated at 0.04lto;o3?, compared to 0.0416lo;oo39 for ^b- 
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I. INTRODUCTION 

The recent discovery of Dark Energy has not only broadened our knowledge of 

the universe, it has brought into sharp relief the degree of our understanding of it. Only a 
small fraction of the mass-energy density of the universe is made up of matter that we have 
characterized; the rest consists of Dark Matter and Dark Energy, both of which have not 
been experimentally detected, and both of whose precise properties are not known. Both 
are needed to explain what is seen on an extremely wide range of length scales. On the 
galactic (~ 100 kpc parsec), galactic cluster (~ 10 Mpc), and supercluster (~ 100 Mpc) 
scales. Dark Matter is used to explain phenomena ranging from the formation of galaxies 
and rotation curves, to the dynamics of galaxies and the formation of galactic clusters and 
superclusters. On the cosmological scale, both Dark Matter and Dark Energy are needed to 
explain the evolution of the universe. 

While the need for Dark Matter is ubiquitous on a wide range of length scales, our 
understanding of how matter determines dynamics on the galactic scale is lacking. Recent 
measurements by WMAP jsl have validated the ACDM model to an unprecedented precision; 
such is not the case on the galactic scale, however. Current understanding of structure 
fonnafo. is based on Q. and both analytical solntions and nun,et,cal si.ulat.ons flB 
|8|, |9|, HOj of galaxy formation have been done since then. These simulations have consistently 
found a density profile that has a cusp-like profile ^, I^, [lO|, instead of the pseudoisothermal 
profile commonly observed. Indeed, De Blok and coworkers |ll| has explicitly shown that the 
density profile from attained through simulation does not fit the density profile observed 
for Low Surface Brightness galaxies; the pseudoisothermal profile is the better fit. 

This is the cuspy-core problem. There have been a number of attempts to solve it 



within ACDM [9|, |10| . with varying degrees of success. While the problem does not exist 
for MOND [l^, there are other hurdles MOND must overcome. Our approach to this 
problem, and to structure formation in general, is more radical; therefore, its consequences 
are correspondingly broader. It is based on the observation that with the discovery of Dark 
Energy, A/j^;, there is a universal length scale, Xde = c/^AdeGY^'^, associated with the 
universe. Extensions of the geodesic equations of motion (GEOM) can now be made that 
will satisfy the equivalence principal, while not introducing an observable fifth force. While 
affecting the motion of massive test particles, photons will still travel along null geodesies. 
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and gravitational lensing is not changed. For a model galaxy, the extend GEOM results in 
a nonlinear evolution equation for the density of the galaxy. This equation is the minimum 
of a functional of the density, which is interpreted as an effective free energy for the system. 
We conjecture that like Landau-Ginzberg theories in condensed matter physics, the system 
prefers to be in a state that minimizes this free energy. Showing that the pseudoisothermal 
profile is preferred over cusp-like profiles reduces to showing that it has a lower free energy. 

Here, phenomena on the galactic scale are inexorably connected to phenomena on the 
cosmological scale, and a cosmological check of our theory is made. The Hubble length 
scale Xh = c/hHo naturally appears in our approach, even though a cosmological model is 
not mentioned either in its construction, or in its analysis. Using the average rotational 
velo city and core sizes of 1393 galaxies obtained through four different sets of observations 
spanning 25 years, we calculate as to be 0.68±o.ii5 in excellent agreement with 
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0.761lQ;Q4g from Q|. We also calculate fiasymp, the fractional density of matter that cannot be 
determined through gravity, to be 0.197±o.oi7, which is nearly equal to the fractional density 
of nonbaryonic matter Qrn — = 0.1961^q o26 We then find the fractional density of 
matter in the universe that can be determined through gravity, i^oyn, to be 0.041!1^q'J]^ which 
is nearly equal to VIb = 0.0416^Q QQ3g. Details of our calculations and theory is in 
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II. EXTENDING THE GEOM AND GALACTIC STRUCTURE 

Any extension of the geodesic action requires a dimensionless, scalar function of some 
property of the spacetime folded in with some physical property of matter. While before 
no such properties existed, with the discovery of Dark Energy there is now A^i^; and these 
extensions can be made. As we work in the nonrelativistic, linearized gravity limit, we 
consider the simplest extension: 

/:Ext = mc(l+D [RcyAnEG] f (^7..^^) ' = mclK[i?cVAz,^G] (^7..^^) ' (1) 

with the constraint v"^ = (? for massive test particles. Here, is a function function 

given below, and R is the Ricci scalar. For massive test particles, the extended GEOM 
is v""^ yV^ = {g'"' - vf^v^/c^) Vy\og9\[A + 8ttT/Aoec^], where -u^ is the four-velocity of 
a test particle, Tf^u is the energy- momentum tensor, T = T^, and we take Ade to be 
the cosmological constant. As the action for gravity+matter is a linear combination of 
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the Hilbert action and the action for matter, any changes to the equation of motion for test 
particles can be accounted for in T^,^, and we still have R = AKdeG / + S'kGT / in Eq. ([1]). 
For massless particles, v'^V y (9^[4 + ^txT j Kj:,ec^]v^) = instead. With the reparametization 
dt the extended GEOM for massless test particles reduces to the GEOM. Our 

extended GEOM does not affect the motion of photons. 

Because the geodesic Lagrangian is extended covariantly, Eq. ([1]) explicitly satisfies the 
strong equivalence principal. For T^,^, we may still take T^u = (p + p/c^)v^Vu — pg^v for an 

Dust 



inviscid fluid with density p and pressure p [16||. While for the GEOM ygeo-Dust _ py^y^ 
for dust, for the extended GEOM the pressure does not vanish [16|; it is a functional of 
p and 9^. Nevertheless, in the nonrelativistic limit p « pc^, and Jj^xt-Dust ^ pv^v^ still 
fiol . Moreover, because v^v^ = (? for the extended GEOM, the first law of thermodynamics 
still holds for the fluid, and the standard thermodynamical analysis of the evolution of the 
universe under the extended GEOM follows much in the same way as before. 

All dynamical effects of extension can be interpreted as the rest energy gained or lost 
by the test particle due to variations in the local curvature. For these effects not to have 
already been seen, !D(4 + SttT/AdeC^) must change very slowly at current experimental 
limits. As such, we take D{x) = x{c(a) I^i^ + s^^"''^)'^ds, where a\ > 1 and x{c(a) is 
set by 2)(0) = 1. This D{x) was chosen for three reasons. First, there is only one free 
parameter, a\, to determine. Second, it ensures that the effects of the additional terms 
in the extended GEOM will not already have been observed; Ade = (7.21+o;84) x 10"^° 
g/cm^, and p ^ A/)^/27r in all current experimental environments so that 2D ^ 0. A 
lower experimental bound of 1.35 for a\ can be found [l^. Third, D'{x) is negative, and 
will contribute an effective repulsive potential to the extended GEOM that mitigates the 
Newtonian 1/r potential. 

While definitive, a first principles calculation of the galactic rotation curves using the 
extended GEOM would be analytically intractable. Instead, we show that given a model, 
stationary galaxy with a specific rotation velocity curve v{r), we can derive the mass density 
profile of the galaxy. We use a spherical model for the galaxy that has three regions. Region 
I = {r I r < Tfj, and p ^ Aeie/^ti}, where th is the galactic core radius. Region II 
= {r I r > Tff, r < rjj, and p ^ Aoe/'^'^} is the region outside the core containing 
stars undergoing rotations with constant rotational velocity; it extends out to rjj, which is 
determined by the theory. A Region III = {r | r > rjj, and p <^ Ae/e/^tx} also appears in 
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the theory. 

As all the stars in the model galaxy undergo circular motion, the acceleration of a star, 
a = X, is a function of is location, x, only. Taking the divergence of the extended GEOM, 

1 [ 2 I + a\ f 8n 



where k\p) = {l + {A + Sirp/ AoEf^""^} /x^Ie^ and /(x) = -V ■ sl/AttG. We do not 
differentiate between baryonic matter and Dark Matter in p. Near the galactic core 1 / /t(p) ~ 
ADi?[ADi?/87rpi^](i+"A)/2^ ^l^gre pn is the core density. Even though Xde = 14010t8?o Mpc, 
because pn ^ Ade/^tt, q;a can be chosen so that l/n{r) is comparable to typical rn- Doing 
so sets a\ ~ 3/2. 

Given a v{r), a(r) can be found and /(r) determined. We idealize the observed velocity 
curves as ■y^'^^^'(r) = Vht/th for r < Th, while ■^'^'''^^(r) = Vh for r > r^, where Vh is the 



observed asymptotic velocity. 



his -yi'i'^ai^r) is more tractable than the pseudoisothermal 



As it has the same limiting forms in both the r <^ r^j 



velocity curve, fP"'^°(r), used in 
and r ^ th limits, ■y"^'^'^^(r) is also an idealization of fP"'^°(r). 

For cusp-like density profiles [10'], it is the density profile that is given. While it is possible 
to integrate the general density profile to find the corresponding curves fcuspl'"), both the 
maximum value of fcusp(^) and the size of the core are different depending on the profile. 
These core sizes would thus have to be scaled appropriately to compare one profile with 
another. Doing so is possible in principle, but would be analytically intractable in practice. 
We instead take /(r) = pu {rn/ry if r < rn, and /(r) = pn (r/f/r)^/3 ii r > rn for the 
density profiles. Here, 7 < 2 and (3 >2 agrees with the parameters for the generic cusp-like 
density profile [3], with the core size set to th- The 7 = 0,/3 = 2 case corresponds to the 
idealized psuodoisothermal profile. 

Since p ^ A£,e/27c in Regions I and II, Eq. ([2]) minimizes 
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which we identify as a free energy functional; here, u = r/x^^'^^DE- For 7 = 0, 
Eq. ([2]) gives p(r) = pn in Region I; the free energy for this solution is ^J^'y=o = 
—Adet^h {-^de/Sttph)"^"^ /6(aA — !)• While for 7 > perturbative solutions can be found, 



all such solutions have a ^JF^ greater than ^T^=o [iG*]. This results because ~ |Vpp > in 
Eq- ©; just as in a Landau-Ginzberg theory, |Vpp only vanishes for the constant density 
solution. 

For Region II, the density, p//, is first found asymptotically in the large r limit. With 
the anzatz /(r) ^ p(r) for large r, Eq. ([2]) reduces to a homogeneous equation [igI with the 
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solution Pasymp(M) = Az)£;S(aA)/87rM i+"A , where S(aA) = [2(1 + 3aA)/(l + aA)T+"A • To 
include the galaxy's structural details, we take pair) = Pasymp('") + pjii^) and to first order 
in pji, 

Pii{r) = pi>symp{r) + -Af3pH y—j + y—j (Ccos cos [z/Q log r/r/z] + C^in sin [z/q log r/r^]) • 

(4) 

where z/q = [2(1 + 3a;A)/(l + ^a)^ — 1/4]"^^^, Ccos and Cgin are determined by boundary 
conditions, and A/j = 1 for /5 = 2,3. The first part, Pasymp('"), of pn{r) corresponds to 
a background density. It is universal, and has the same form irrespective of the detailed 
structure of the galaxy. The second part, P//(r), gives the structural details. 

The free energy, ^^JF, for Region II separates into the sum of three parts. The first part 
depends only on Pasymp; it is positive, and is independent of (3. The second part is 

-o'l d'^fi-) ( ^^y^+ . 'ZnL. [ «V}.(«)Vpasymp-rfS, (5) 
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{X^'^^DeY Jdj, ' ' VSTTPasymp; Afli^S^d+^A) 

where Djj is Region II. It is negative because the minimum p must be positive. Indeed, we 
find that ^^^asymp-/3 ~ -(r/f/r,,)/^ for P < 5/2; ^^^asymp-/3 ~ -{rH/r-nf/^ for 5/2 < /5 < 
5 - 2/(1 + «a); and ^^^asymp-/3 ~ ±(rH/x^/'AD,,)5-2/(i+"A) for 5 - 2/(1 + c^a) < /5. Clearly, 
free energy is lowest for (3 = 2. The third part depends on {p}j{r))'^, and is negligibly small. 
The total free energy in this region is thus smaller for (3 = 2 than for (3 > 2. Combined with 
the calculation for ^JF, we conclude that the pseudoisothermal density profile has the lowest 
free energy, and is the preferred state of the system. We thus take 7 = and (3 = 2 in the 
following. 

In Region III, p < ADE/2n, and K^{p) ^ (1 + 4i+"^)/xAf)^; Eq. reduces to the 
undriven, modified Bessel equation. As such, the density vanishes exponentially fast in this 
region on the scale l//t(r). This sets rjj = [x/(l + 4:^~^"^)Y^'^^de- 

The extended GEOM can be written as x = — V2J. The dynamics of test particles is 
governed by an effective potential 2J(x) = $(x) + c^log(!H[4 + STrpA^^;]), and not by the 
gravitational potential $(x). For $(r) in Region I, we obtain the Newtonian gravity result 
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$(r) = v'fjr'^/2rfj + constant. In Region II, $(r) is dominated by four terms. The first is 
the usual 1/r term. The second is a log{r/rH) term due to f{r). This term is long ranged, 
and in addition to galactic rotation curves, could explain the interaction observed between 
galaxies and galactic clusters. The third is a p]j{r)r'^ term, and contains terms ~ 
The fourth term is a r^"'^/^^+"*) term due to Pasymp, and is proportional to c^. 

This last term grows as r^/^ for oa = 3/2, and would dominate the motion of test particles 
in the galaxy if the extended GEOM depended on $(x) instead of 23(x). We instead find 
that 5J(x) ^ $(x) - [u'^xc^i^ + «A)^/4aA(l + 3a;A)] (AD£;/87r)(pasymp - ttAP//)- The last two 
terms in this expression cancel both the p\i{r)r'^ and the r^°^/*^^+"^^ terms in $(r); the 
resultant 23(r) increases as logr/ru, agreeing with observation. 

The r'^^^/^'^+°''^) term in $(x) comes from the background density Pasymp- Thus, a good 
fraction of the mass in the observable galaxy does not contribute to the motion of test 
particles in the galaxy. It is rather the near-core density p}/(r) that contributes to 23(x). 
As inferring the mass of structures through observations of the dynamics under gravity of 
their constituents is one of the main ways of estimating mass, the motion of stars in galaxies 
can only be used to estimate p\j] the matter in pasjmpij) is present, but cannot be "seen" 
in this way. Moreover, as Pasympl'") ^ P//(^) when r ^ ru^, the majority of the mass in the 
universe cannot he seen using these methods. 

III. A COSMOLOGICAL CHECK 

We have extrapolated our results for a single galaxy to the cosmological scale. This is 
possible because recent measurements from WMAP, the Supernova Legacy Survey, and the 
HST key project show that the universe is essentially flat; h = 0.732^q'q32 and of the age 
of the universe to = 13.73^g Gyr were determined using this assumption. The largest 
distance between galaxies is thus do = K{Q)Xh, where = 1.03±o.o5- 

Next, the density of matter of our model galaxy dies off exponentially fast at rjf, the 
extent of matter in the galaxy is fundamentally limited to 2r//. This size does not depend on 
the detailed structure of the galaxy; it is inherent to the theory. Given a I^a = 0.716-|-o.o55) we 
can express r// = 

[87rx/3fiA(l + 4i+"A)]V2^^ [Igj as well Q, and numerically ru = 0.52Xh 
for ftA = 3/2. Although c^a was set to 3/2 based on analysis at the galactic scale, p(r) 
naturally cuts off at Ai//2. 
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To accomplish the extrapolation, we consider our model galaxy to be the representative 
galaxy for the observed universe. This representative galaxy could, in principal, be found by 
sectioning the observed universe into three-dimensional, non-overlapping cells of different 
sizes centered on each galaxy. By surveying these cells, a representative galaxy, with an 
average and r|^, can be found, and used as inputs for the model galaxy. Even though 



such a survey has not yet been done, a large repository of galactic rotation curves and core 
[isl is present in the literature. Taken as a whole, these 1393 galaxies are 
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reasonably random, and are likely representative of the observed universe at large. 

While we were able to estimate of q;a = 3/2 by looking at the galactic structure, the 
accuracy of this estimate is unknown; comparison with experiment is not possible. We 
instead require that rjj = M{Q)Xh/2, which in turn gives a\ as the solution of ^{QY{1 + 



41+aA) ^ 327rx(«A)/3fiA; this sets 
A calculation of as has been done 



•51±o.ii- 

16| using Eq. (jlj). The resultant (t| is dominated by 
two terms. The first is due to the background density Pasymp- It depends only on a^, and 
contributes a set amount of 0.141 to cr|. The second is the larger one, and is due primarily 
to the 1/r^ term in Eq. (jl]). It depends explicitly on the rotation curves through the term 
K/c)4(8/i-iMpc/r|,). 

Although there have been a many studies of galactic rotation curves in the literature, 
both vh and r^ are needed here. This requires fitting the observed velocity curve to some 



model. To our knowledge, both values are available from four places in the literature: 



111]; the CP data set [14]; the Mathewson et. al. data set 
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'he 



de Blok et. al. data set 

analysed in jl^; and the Rubin et. al. data set [13]. Except the last set, the observed velocity 
curves is fitted to either t>P"^'^°(r), or to a functionally similar velocity curve [l^. The last set 
gives only the galactic rotation curves, and they have been fitted to v^~^^°{r) in |l6|. While 
the URC of [18] has a constant asymptotic velocity, it has a r°'^^ behavior for r small. This 
behavior is different from ti'^^^i^ and was not considered here flG ]. 

While Vh is easily identified for all four data sets, determining rn is more complicated; 
this is determination is done in 1^. The resultant values are used to obtain and for 
each set, which are then used to calculate the Ug and Acrg for it. Results of these calculations 
are in Table [B Four of the five data sets give a dg that agrees with the WMAP value at 
the 95% CL. The Rubin et. al. set does not, but it is known that these galaxies were not 
randomly selected isl ]. 
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Data Set 




11 


11 


Ar* 

il 


^8 


Acts 


t-test 


deBlok et. al. (53) 


119.0 


6.8 


3.62 


0.33 


0.613 


0.097 


1.36 


CF (348) 


179.1 


2.9 


7.43 


0.35 


0.84 


0.18 


0.43 


Mathewson et. al. (935) 


169.5 


1.9 


15.19 


0.42 


0.625 


0.089 


1.34 


Rubin et. al. (57) 


223.3 


7.6 


1.24 


0.14 


2.79 


0.82 


2.46 


Combined (1393) 


172.1 


1.6 


11.82 


0.30 


0.68 


0.11 


0.70 



TABLE I: The (km/s), (kps), and resultant as, Acrg, and t-test comparison with the WMAP 
value of as ■ 

We have estimated fiasymp by averaging Pasymp('^) over a sphere of radius rjj, and found 
f^asymp = 0.197±o.oi7- In Calculating this average, we assumed that there is only a single 
galaxy within the sphere, however. While this is a gross under counting of the number 
of galaxies in the universe, Pasymp is an asymptotic solution, and p]j rapidly with r. 
Additional galaxies may change the form of Pasymp, but these changes are expected to be 
equally short ranged; we expect that our calculation is an adequate estimate of fiasymp- 
Such is not the case for flByn, however. Direct calculation of fioyn would require knowing 
both the detailed structure of galaxies, and the distribution of galaxies in the universe. 
Instead, we note that ilm = ^asymp + ^Dyn, and using — 0.238_Qg2g from WMAP, find 

^Dyn = 0.0411q;q3°. 



IV. CONCLUDING REMARKS 



Given how sensitive as is to v^, r*^, and oa, that our predicted values of as is within 
experimental error of the WMAP value is surprising. Even in the absence of a direct exper- 
imental search for oa, this agreement provides a compelling argument for the validity of our 
extension of the GEOM. It also supports our free energy conjecture; our calculation of as 
would be very different if /5 = 3, say, was used instead of /5 = 2. With = 1.51 so close to 
the experimental lower bound for a\ of 1.35, direct measurement of may also be possible 
in the near future. 

Interestingly, fi„ - I^b = 0.196l|]:|J| is nearly equal to fiasymp in value. Correspondingly, 
fis jsl is nearly equal to fioyn- It would be tempting to identify fiasymp with fi^ — ^b, espe- 
cially since matter in Pasymp('") is not "visible" to inferred- mass measurements. That Qnyn 
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would then be identified with Jl^ is consistent with the fact that most of the mass inferred 
through gravitational dynamics are indeed made up of baryons. We did not differentiate 
between normal and dark matter in our theory, however. Without a specific mechanism 
funncling nonbaryonic matter into Pasymp and baryonic matter into p — Pasymp, we cannot at 
this point rule out the possibility that — VLb = f^asymp and VLb ^ i^Dyn is a numerical 
accident. 
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